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We study next-to-leading order contributions to the soft static fermion dispersion relation in hot 
QED. We derive an expression for the complete next-to-leading order contribution to the retarded 
fermion self-energy. The real and imaginary parts of this expression give the next-to-leading order 
contributions to the mass and damping rate of the fermionic quasi-particle. Many of the terms 
that are expected to contribute according to the traditional power counting argument are actually 
subleading. We explain why the power counting method over estimates the contribution from these 
terms. For the electron damping rate in QED we obtain: 'jqed = ^^{2.70). We check our 
method by calculating the next-to-leading order contribution to the damping rate for the case of 
QCD with two flavours and three coulours. Our result agrees with the result obtained previously 
in the literature. The numerical evaluation of the nlo contribution to the mass is left to a future 
publication. 

PACS numbers: ll.lO.Wx, 11.15.-q 



I. INTRODUCTION 



It is well known that the behaviour of elementary particles becomes modfied when the particle propagates in a 
medium. The particles become "dressed" by their interaction with the medium, and one speaks of collective modes, 
or quasi-particles. One studies these collective modes by looking at the corresponding thermal propagators. The 
behaviour of the quasi-particles is deduced from the analytic structure of the propagator. In Q it was shown from 
general principles that the singularity structure of certain components of gauge and matter propagators are gauge- 
independent, when all contributions of a given order are systematically taken into account. 

However, the actual calculation of dispersion relations for soft quantities at next-to-leading order is notoriously 
difficult. The original paper by Braaten and Pisarski Q identified three potential contributions at next-to- leading 
order. They are: (1) corrections to the result for the 1-loop diagram obtained by expanding to next-to-leading order 
in the ratio of the external momentum to the loop momentum; (2) contributions to the 2-loop diagrams from the 
region of the phase space that corresponds to both loops hard; and (3) the contribution to 1-loop diagrams with soft 
loop momentum, and all propagators and vertices replaced with htl effective ones. We will refer to these three types 
of contributions as (1-terms), (2-terms), and (3-terms). We note that the power counting arguments of Braaten and 
Pisarski refer to the maximum possible contribution from each type of term. The actual contribution may be lower 
order for kinematical reasons, or because of some cancellation between different integrals. 

Only a few calculations of next-to-leading order quantities have been done. The gluon damping rate was calculated 
by Braaten and Pisarski Q. The damping rate of a soft static quark was calculated by Kobes, Kunstatter and Mak 
[j| , and by Braaten and Pisarski Q . All three of these calculations involve only the imaginary part of the appropriate 
2-point vertex function. The real part is generally more difficult to calculate. The gluonic plasma frequency in the 
long wavelength Hmit has been calculated by Schulz 0|, and the Debye mass has been studied in p]. 

In this paper we consider the next-to-leading order corrections to the dispersion relation for soft static fermions 
Qfj. = {<la ~ eT, g — 0). For static fermions, the collective modes appear mathematically as poles of the propagator. 
The mass and damping rate of the quasi-particle are obtained from the solution of the equation: 

det iQ - SretiQ)) =0. (1) 

qo=M—ii 

The leading contributions are obtained from the hard thermal loop (htl) results: 



ruf = ^e^rys ; t^") = . (2) 

The goal is to find an expression for the next-to- leading order corrections to this result. 

One of the main results of this paper is a compact expression for the next-to-leading order contributions to the 
retarded fermion self-energy in QED. This result can be divided into real and imaginary parts (Eqns. II33|I and (EUl) 
which give, respectively, the next-to-leading order corrections to the mass and damping rate in The expression 
for the imaginary part can be checked by converting to the corresponding QCD expression (obtained previously in 
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B'sl) by calculating the appropriate group factors. The expression for the real part has not appeared previously in 
the literature. 

The imaginary part of the next-to-leading (nlo) self-energy gives the nlo contribution to the damping rate 



^Tr(7oImS(go,0)) . (3) 

qo=mf 
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To check our method, we have calculated the nlo quark damping rate for QCD with 2 flavours and 3 colours and 
obtained a result that agrees with that of (Note that these authors use a different convention for the definition 

of the damping rate: they define the damping rate with an extra factor of 1/2 relative to the definition used here). 
For QCD with 2 fiavours and 3 colours we obtain: 

{Nf = 2, N, = 3} ^ j=^^^i2.81); Cj. = 4/3 . (4) 

47r 

The corresponding result for QED has not been previously calculated. The result is: 

IQED = ^{2.70) . (5) 

The calculation of the real part of the fermion self-energy is more complicated for several reasons which we discuss 
below. 

(1) There are pure real tadpole type contributions that can be dropped in the calculation of the imaginary part. 

(2) The leading order contribution to the imaginary part of the self-energy is zero (because of kinematic constraints), 
and thus the complete evaluation of the integral gives a contribution of e^T times some number of order one, which 
is the full nlo contribution. For the real part, the numerical result for the integral will contain the leading order term 
(which is of order eT), in addition to the nlo term. The leading order term must be subtracted as a counter-term to 
obtain the numerical result for the nlo contribution. 

(3) For the real part of the electron self-energy, the (1-terms) have been calculated in Q and shown to be of 
order e'^Tln(l/e), and not of order e^T. In addition, we have shown by explicit calculation that the contribution 
from (2-terms) is also of order e^Tln(l/e). These results appear to indicate that the power counting argument over 
estimates the contribution from (1-terms) and (2-terms) terms. A similar result was obtained by Schulz in f^: he 
found that the contribution to the gluonic plasma frequency from 2-loop diagrams with hard internal momentum is 
smaller than predicted by power counting. It is not clear if (2-terms) terms are always over estimated by the power 
counting argument. This issue is discussed in Appendix E. 

The ln(l/e) part of the contribution discussed above comes from an integral of the form 

dpp {mip) + nf{p)) Prin^-^— 2 (6) 

P ~ Qo 

where the notation "Prin" indicates a principle part and the factors rihip) and nf(j)) are bose-einstein and fermi-dirac 
distribution functions. The integral is cut off at large T by the distribution functions and regulated at small p by the 
soft cutoff go- The point is that both of the distribution functions contribute to the result. This is not the case for 
the nlo contribution to the imaginary part which is dominated by the soft part of the momentum integral and the 
infra-red part of the bose-einstein distribution function. The result of the calculation in ^ |^ justifies dropping all 
terms proportional to the fermi-dirac distribution function in the integrand for the imaginary part. In the expression 
for the real part we keep all terms containing either bose-eintein or fermi-dirac distributions. 

(4) The numerical calculation of the real part of the self-energy is also more complicated than the corresponding 
integral for the imaginary part. The integrand for the real part contains a cut that gives large positive and negative 
contributions as it is approached from either side. These limits must be handled carefully to extract the surviving 
finite contribution. 

A compact analytic result for the real part of the fermion self-energy at nlo is given in Eqn. II34II . The result of 
the numerical calculation of this integral will be presented in a future publication. 

This paper is organized as follows. In section ^ we define some notation and discuss the diagrams that contribute 
to the fermion self-energy at nlo. In section IIIII we give some information about the technical aspects of the finite 
temperature calculation. In section IIVI we define our notation for the propagators and vertices. In section El we 
describe the structure of the calculation and give the analytic results for the integrands of the retarded self-energy, 
and its real and imaginary parts. In section IVll we give the numerical result for the imaginary part. In section IVlTI 
we present our conclusions. In Appendix A we explain the notation we use to assign indices to the components of 
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the Keldysh n-point functions and vertices. In Appendix B we give the full expressions for the discontinuities of the 
propagators. In Appendix C we Hst the ward identities and kms conditions that are needed to do the calculation. In 
Appendix D we give some explicit results for the vertex functions that cannot be rewritten using ward ientities. In 
Appendix E we discuss power counting arguments. 



II. PRELIMINARIES 



We use capital letters to denote 4-momenta: K — {ko,k). We take the external momentum to be Q = {qQ,0) and 

2 
/ 



write R = P + Q so that we have r= p. The fermion and photon thermal masses are: rrii = and uiq — 



We use: 

{7m' 7i^} = 2g^i> ; 9t^iy = diag (1, -1, -1, -1) . 
The thermal distribution functions are defined as: 

""'^P^^^^Firj'^ "/(p) = ;^;Vt' ^sb) = l + 2nfc(p); NF{p)^l-2nf{p). (7) 

In this paper we are only interested in thermal effects and consequently we ignore zero temperature pieces of the 
self-energy. The fermion self-energy is decomposed in the usual way: 

S = 7°I]W+7-gS«. (8) 

Since we have taken q — the only non-zero component is I]'-*'-' — Tr(7°S)/4. From now on, we suppress the 
superscript '(0)' to simplify the notation. We use throughout the conventions of but we work in real time. 

The diagrams to be calculated are shown in Fig. 1, where dotted fines indicate photons and sofid lines are fermions. 
The notation is defined in Fig. 2 where P^in indicates the momentum of an incoming fermion, P^out is the momentum 
of an outgoing fermion, and P!^^^ is the momentum of an incoming photon. The vertices are given by: 

r^C^Vin, -fVout) = 7m + rhtlp(P0iii, P^out) (9) 
-^/ii^(PV'in7 Pjin^ ^'yin^ P^out) — -^'-^htl^ix (P^jn , P^ij^i ^^ini P^ont) 

The htl vertex functions can be written as the contraction of a gamma matrix and a factor that does not contain 
dirac structure. We introduce the notation: 

Ehtl := 7^EhtlA (10) 



rhti^ := 7^rhti^ 

Afhtl^^ := 7^Afhtl^^A 



In order to keep the notation simple, we will not add a tilde to the vertices on the rhs of these expressions. The extra 
lorentz index indicates that the gamma matrix has been factor out. 



P/ \ P ^ 



Q R=p+Q Q Q Q 



FIG. 1: The diagrams to calculate for a soft electron. 
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i e M (Pvm,??^ ,PYin ,Pvou) 



FIG. 2: Definitions of conventions for propagators and vertices. 

All of the contributions identified by the power counting arguments of Braaten and Pisarski as contributing at nlo, 
as well as the leading-order contributions, are formally included in Fig. 1 if the integral over the internal momentum 
variable is not restricted to any particular region of the phase space, but runs over the full range of the momentum 
variable. 



III. REAL TIME FINITE TEMPERATURE FIELD THEORY IN THE KELDYSH REPRESENTATION 

We will do the calculation in the Keldysh representation of the real time formulation of finite temperature field 
theory. The basic structure of the method is described below for a scalar theory. Details are given in In the 

closed time path (CTP) formulation of finite temperature field theory, the path integral is defined on a contour in 
complex time that has two branches. The '1' branch runs from minus infinity to positive infinity just above the real 
axis, and the '2' branch runs back from positive infinity to negative infinity just below the real axis. 



Ci 






Bs» 1 


< 





FIG. 3: The CTP contour in the complex time 
Green functions are defined on this contour: 

iDc{x,y) = {fc(j){x)(t){y)) 

where indicates time ordering along the contour. In the 1-2 representation, a propagator with real time arguments 
is written as a 2x2 matrix: 

iDnix,y) = {T(b{x)(biy)) 
iDi2{x,y) = {(j){y)(j){x)) 
iD2i{x,y) = {cl){x)(j){y)) 
iD22{x,y) = {T*<j>{x)(j>iy)} 
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where the operators T and T* indicate time ordering and anti-time ordering, respectively. Note that there are at 
most three independent components because of the constraint 

Dii~Di2- 021+0-22^0. (12) 

In general, in the 1-2 representation, an n-point function is a tensor with 2" components, corresponding to the two 
possible values of each index. Only 2" — 1 of these components are independent because of a constraint of the same 
form as the one given above. Truncating external legs gives the corresponding vertex functions. There are (2^ — 1 = 3) 
self-energy functions, (2'^ — 1 = 7) 3-point vertices and (2'* — 1 = 15) 4-point vertices. 

In the Keldysh representation, we use linear combinations of components in the 1-2 representation that have a more 
direct physical interpretation. For example, the three components of the propagator in the Keldysh representation 
correspond to the retarded, advanced, and symmetric propagators. For any Keldysh vertex function, each of the 
2" — 1 components is denoted by a single index which takes values {2, 3, 4, • • • 2"}. The definitions we use for assigning 
indices to the Keldysh functions are given in Appendix A. For clarity, where possible, we will write propagators and 
self energies as retarded, advanced, or symmetric functions instead of using the corresponding numerical index. 

We give explicitly below the expressions for the three propagator components, and the retarded and advanced self 
energies. The retarded, advanced and symmetric propagators are given by: 

Oadv{K) = On - O21 = ^ 



X2 - zSign(fco)e 

0,y^{K) - Dn + D22 = NB{kQ){Oret{K) - Oadv{K)) = -2^ziVB(A:o)Sign(fco) 5{K'') . 

The retarded and advanced self energies (which give the real and imaginary parts: ReSret = l/2(Sret + 
^adv) ; ImEret = l/(2i)(Sret - ^adv)) are: 

Sret = -f- S12 ; Sadu = '^11 + 5^21 • (14) 

The calculation of the diagrams in Fig. 1 is performed as follows. Each line carries an index on each end which can 
take values 1 or 2. One sums over all internal indices and takes the combinations of external indices that corresponds 
to the component of the self-energy that is wanted. The result is rotated to the Keldysh basis. The technique is 
described in [lH . 

IV. HTL PROPAGATORS AND VERTICES 

The basic structure of QED at finite temperature is the same as for scalar theory, but we have to include the dirac 
and lorentz structure of the propagators and vertices. 

A. Propagators 

In Eqns. IjlRII . Illfill and ()17ll . we simplify the notation by supressing the subscripts ret, adv. For photons the 
symmetric propagator is given by 0'^y^{P) — NB{po){0'^^f{P) — D'^'^^{P)) and for fermions we have Ssym{R) — 

NF{ro){Sret{R)- Sadv{R))- 

The fermion propagator is written 

S{R) = - Shtl(i?))-^ (15) 
The inverse propagators are inverted by writing 

S{R) = i(5+(i?)(7o - f) + ^5_(i?)(7o + f)) (16) 



S+{R) 



2'^' " 2 

2r2 



(mj + r{r- m)) + In (r - ro) mj 



2r 



9r 

S-iR) 



2r (mj + r{r + tq)) - In (^^'^ mj (r + ro) 
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We use the covariant gauge and write the photon propagator in terms of transverse and longtitudinal projections 
(recall that p = r): 



D,AP) - P/.MP) + Pf:.^DL{P) (17) 

p2 



Dt{P) = ^^^^, r; Dl{P) 



G(po,?-) = ^|l JPrnQ + mQ; Qo(Po,r) 

F{po, r) = — 2 2™G 1 



^ In ' ^° ^ 



2 VPo - 



Pl^,^^.<rQl,L(p._e^) (p._£^ 



P^ qi 90 A 90 

The expressions for the projection operators have been obtained by writing the vector that defines the rest frame of 
the medium as — {1,0,0,0} — Q^/q^. 

The discontinuity of a scalar propagator is defined as: 

p{K) = -I d{K) = -l{Dret{K) - Dadv{K)) (18) 

The discontinuities of the fermion and photon propagators are obtained from: 

d+ (R) = S+, (i?) - 5+ , (R) ; d_ (R) = a,- , (i?) - S-,^ {R) (19) 
driP) = Dl,{P) - D^^,{P) ; d^iP) = D^P) ^ D^ad.{P)) ■ 

Each discontinuity has a pole and a cut contribution and has the general form: 

da,(K)^-2TTi Za,(k)nS{ko~nuj^{k))-2Tri(3a,{ko,k); x = {T,L,+,-} (20) 

ri=±l 

The functions ujt{p), ^l{p), ^+{p) and i^--{p) are the roots of the denominators of the corresponding propagators. 
The expressions for these functions, and for the residues Zi and cut functions Pi are given in Appendix B. 

V. CALCULATION OF THE HTL CORRECTED 1-LOOP DIAGRAMS 

A. structure of the calculation 

The first step in the calculation of the diagrams in Fig. 1 is to sum over indices in the 1-2 representation and rotate 
the result to the Keldysh basis. We define the operator 

which will be factored out of all expressions. We use the techniques of to obtain: 

^retiQ) (22) 
diagram 1 



= TV Tr 



^o^sym(p)p. [ 2 ] (Q, R)Sret {R)r^ [ 2 ] (i?, Q) + {P)T'' [ 2 ] (Q, R)S,^^{R)T^ [ 3 ] (i?, Q) 



Bf: {P)V- [ 2 ] (Q, R)Sr,t {n)V^ [ 4 ] (i?, Q) + i^^t* (^^)r'^ [ 4 ] (Q, R)S,,t {R)V^ [ 2 ] (i?, Q) 
{PW [ 6 ] (Q, R)Saa. (i?)r^ [ 3 ] (i?, Q) 
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2 



(23) 

diagram 2 

7° (i?^r (^)A^^''[2](Q, P, ~P, Q) + i^;:*(P)M'^''[4](g, P, -P, Q) + Z?^f (P)Af^''[6](Q, P, -P, Q)) 



The indices in square brackets refer to the Keldysh components of the vertices. The notation is defined in Appendix 
A. 

The basic strategy of our calculation is to rewrite the htl vertices in terms of the self energies, and rearrange the 
result to obtain the simplest possible form. In order to do this, we first note that there are four basic structures in 
the transverse and longtitudinal projection operators ifTTjl : 

P.P.; {^p.,^p4; 5.- (24) 

The first two forms give a zero contribution to the final result. This is a consequence of gauge invariance and can be 
derived using the ward identities. One finds that the gauge dependent contribution to the nlo fermion self-energy is 
proportional to an integral times the square of the inverse propagator S'^^^^(Q) — Q — Shtl((5), which vanishes on the 
mass shell. In it was pointed out that a straightforward evaluation of the integral produces mass-shell singularities 
that cancel the contributions from the two inverse propjigators and give a finite gauge dependent contribution to the 
damping rate. This problem was resolved by Rebhan who showed that the integral must be regulated before the 
mass shell is approached. Using this procedure one finds that the position of the pole is gauge independent, and the 
gauge dependence occurs only in the unphysical residue. 

The third and fourth forms in Ij24|l give non-zero contributions. All terms that come from the third form (Q.Qi//(?o) 
can be rewritten in terms of the fermion self-energy by using the htl ward identities. Terms proportional to g^i, are 
more difficult to handle. We use the kms conditions for the 3-point functions, and some equations obtained by taking 
the complex conjugate and interchanging the arguments of the fermion legs, to obtain the simplest possible expression. 
A complete list of the ward identities and kms conditions for 3- and 4-point functions are found in [^. In Appendix 
C we have collected the ones that we use in this calculation. 

In order to express the results in a compact form, we rewrite the photon propagator in itTzjl keeping only the terms 
that will give non-zero contributions: 

D',AP) = 9,.Dt{P) + ^4^{bT{P) + Dl{P)) (25) 

% 

Dt{P) = -^Dt{P) ; Dl{P) = -Dl{P) ■ 

We will use the same notation for the discontinuities as before: dx/LiP) — -C^t/l(P) ^ ^t/K^)- Using the fact that 
the htl 4-point vertex satisfies Afhtl^'' = 0, we consider the structure of the terms that result when we use lf25|l in 
ll22)l and lf23jl . There are three types of terms: 

(1) contributions from the first term in Ij25|l 

(2) contributions from the second and third terms in ll25ll which can be divided into two types: 
(2a) terms which do not contain fermion propagators 

(2b) terms that do contain fermion propagators 

All contributions to Y^ret{Q) have one of the forms below: 

(l) = 7Vr2dT(P)A^s(po)Tr(7V^,tt(P)7^)rrM[2](g,i?)r/[2](g,P) (26) 

(2a) ^-\n r^dT/L{p)M-f°j^) (2Qr + Rr) Nb (pq) (27) 
% 

(26) = l^N {Qr + Rr) Tr [7%^ (dT,L{P)StAR)NB (po) + D!¥'/l{P)d^{R)NF (ro)) 7^1 (Qa + Rx) (28) 

L \ /J 

Using Eqns. (EJ and itTo)! we have written: 

rr^[i]{Q,R)^gr^. + rht\r^[i]{Q,R); z = 2 or 5 . (29) 
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There are two contributions from terms of type-(l) which come from the (+/ -) modes of the electron. There are two 
contributions from terms of type- (2a) which come from transverse and longtitudinal photon modes. There are four 
type- (2b) contributions which come from the (transverse/longtitudinal) x (+/ -) modes of the photon and electron. 

B. Longitudinal results 

Longtitudinal modes get contributions from terms of type-(2a) and type-(2b) only. For future convenience, we 
also divide results into contributions that depend on boson distribution functions and contributions that depend on 
fermion distribution functions. The contributions to Eret(Q) from longtitudinal photon modes are: 

Type2a(Long, A^b) = —Nr^dL{P)NB (po) (30) 
Type2a(Long, A^i?) = 

Type2b (Long, +, TVs) = - ^ iVr^ (-r + qo + ro f dL{P)S+,t{R)NB (po) 

% 

Type2b (Long, -, TVs) = iVr^ (r + go + r^f dL{P)S;,t{R)NB (po) 

% 

Type2b (Long, +, TVj.) = - ^ TV (-r -f go + ^o)' al (P)d+ {R)Nf (ro) 

% 

Type2b (Long, - , iVj.) = - A TV (r -f go + ro)' aL(P)rf- {R)Nf (ro) 

% 

C. Transverse results 

Transverse modes also contain contributions from terms of type-(l). These terms can be further simplified by using 
explicit results for the htl 3-point vertex functions, which have a particularly simple form when one of the fermions is 
not moving. The expressions that we need are given in Appendix D. Using these results, all components of the vertices 
in ll26ll can be written as simple functions of the htl self-energy. These self-energies also appear in the denominators 
of the htl fermion propagators. The general strategy is to rearrange terms in the numerators so that we can cancel as 
many terms as possible with the corresponding terms in the denominators. Significant simplifications occur when we 
combine terms and use the mass shell condition go — mj. We define the notation: 

Ehtlo(s)(i?) = I]htlS"*(i?) + Shtlg'*''(i?) ; Shtlo(d)(i?) = EhtlS'*(i?) - J^htlf^iR) . (31) 
The contributions to Sret(Q) from transverse photon modes are: 

Type2a (Trans, iVs) = dT{P)NB (po) {~3roR^ + imjro + St^qo + 6P^qo - i?'EhtlS'=*(i?)) 

% 
2 - 

Type2a (Trans, Np) ^ — N {r^ - r^) aT{P)NF (ro) Shtlo(d)(i?) (32) 
% 

Type2b Cftans, +,Nb) = -^N {po - rf (r +po + 2go)' dT{P)NB (pa) S+.tiR) 

% 

Type2b (TVans, - , TVs) = -\ N {r + pof {-r + pa + 2qof dT{P)NB (po) S;,t{R) 

% 

Type2b Cftans, +,Nf)^-\n {po - rf {r + po + 2qof aT{P)d+{R)NF (ro) 

% 

Type2b (TVans, TV^) = -\n {r+ pof (-r + po + 2qof aT{P)d- {R)Nf (ro) 

% 
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D. Imaginary parts 

The damping rate is determined from the imaginary parts of the expressions in 113011 and Ij82ll . The nlo contribution 
comes only from terms that depend on the bose-einstein distribution function. The contributions to Iml^retiQ) are: 

' hr + qo + rof dL{P)d+{R)NB (po) (33) 



Im(L,+,26) 


I 


Im(L,-,26) 


i 


Im(r, 2a) = 


% 


Im(r, +, 26) 


i 


Im(T,-,26) 


i 



N {po - r f (r +P0 + 2qof d+{R)dT{P)NB (po) 



N ir+ pof i-r +P0 + 2qof d-{R)dT{P)NB (po) 



E. Real parts 

As discussed in the introduction, the real parts could receive contributions from terms with both bose-einstein and 
fermi-dirac distribution functions. In addition, there are tadpole-type terms that contribute to the real part. We use 
the notation "Prin" to indicate a principle part. The contributions to ReSret(Q) are: 

Re(L, Nb, 2a) ^ — N r^'dL{p)NB (po) (34) 
Qo 

Re(L, +,NB,2b)^-^ N {-r + qo + ro f dL{P)NB {po) Prin+(i?) 
% 

Re(i, TVs, 2b)^-^ Nr^r + qo+ rof dL{P)NB (po) Prin_(i?) 
1o 

Re(L, +, Nf, 26) = -4 TV (-r + qo + rof d+{R)NF (ro) PrinL(P) 
Qo 

Re(L, -, TVf, 26) = TV r2 (r + go + ro)' d^iR)NF (ro) PrinL(P) 
% 

Re(r, TVs, 2a) = -\ N dT{p)NB [po) (2 {~iroR^ + Amjro + Sr^qo + &P'^qo) - R^TM\o{s){R)) 
% 

Re(r, Nf, 2a) = ^ N R^Np (ro) PrinT(/')Shtlo(d)(i?) 
% 

Re(r, +, iV;^, 26) - TV {po - rf {r + po + 2qof dT{P)NB (po) Prin+(i?) 
% 

Re(r, -, TVs, 26) = -1 TV (r + pof {-r + po + 2qof dT{P)NB (po) Prin_(i?) 
% 

Re(r, +, iV;^, 26) = -\n {po - rf (r + po + 2qof d+{R)NF (ro) PrinT(P) 

Re(r, -, Nf, 2b) ^ -\ N {r + pof (-r + po + 2qof d_ {R)Nf (ro) PrinT(P) 
% 

VI. NUMERICAL EVALUATION OF THE IMAGINARY PART 

The damping rate can be calculated from the imaginary part of the self-energy using © . Each of the four type- 
(2b) terms in Il38|l contains two discontinuities in the form of a product dj'/ ]^{P)d±{R). Each of these discontinuities 
contains an cut piece and a pole piece, as given in Eqns. ijBlll . ljB2ll . IIB3ll and l|B4|l in Appendix B. The type-(2a) 
term in contains one discontinuity of the form driP) which has a cut piece and a pole piece. Thus, the total 
number of terms to be calculated is: (4 x 4) + (1 x 2) = 18. 
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We scale all dimensionful parameters with the leading order fermion mass: rp — > rufro; r mfr\ uiq = xmj. The 
result for the quark damping rate for 2 flavours and 3 colours is: 

^=^^^^^^^; CF = {N^-l)/i2N,): {7Vc-3, Nj = 2} =^ x = 4 (35) 

7qcd(x = 4) = ^^(2.81) 

This result agrees with that of 0, 0| . 

For the electron damping rate in a QED plasma we find: 

X = ml /ml = 4/3 ; ^ -/qed = ^(2.70) (36) 

477 



VII. CONCLUSIONS 



We have studied next-to-leading order contributions to the soft static fermion dispersion relation in hot QED. We 
have derived an expression for the complete next-to-leading order contribution to the retarded fermion self-energy. The 
real and imaginary parts of this expression give the next-to-leading order contributions to the mass and damping rate 
of the fermionic quasi-particle. Many of the terms that are expected to contribute, according to the traditional power 
counting arguments, are actually subleading. We have explained why the power counting method over estimates the 
contribution from sub-leading pieces of the 1-loop diagram, and each of the 2-loop diagrams with both loop momenta 
hard. ^ 

For the electron damping rate in QED we obtain: jqed = %|-(2-70). We have checked our method by calculating 
the next-to-leading order contribution to the damping rate for the case of QCD with two flavours and three coulours. 
Our result agrees with the result obtained previously in the literature. The numerical evaluation of the nlo contribution 
to the mass is left to a future publication. 



APPENDIX A: DEFINITIONS OF THE KELDYSH VERTICES 



As explained in section Hill green functions and vertices in the Keldysh representation are Hnear combinations of 
the green functions and vertices in the 1-2 representation. An n-point function in the 1-2 basis is written: Cfc^, 62, f)3, 
where the indices bi and take the values 1 or 2. Vertices are obtained by truncating external legs: 

Keldysh indices are written and are assigned the values a — 1 :— r and a ^ 2 := a. An n-point function in the 
Keldysh basis is written: Ga^, 02,03, ---an- Vertices are obtained by truncating external legs: 

To rotate from the 1-2 representation to the Keldysh representation we use the transformation matrix: 

The n-point function in the Keldysh representation is given by: 

Go,...a„ = 2t-i (7„^ • • • [/„„''"Gb,...fc„ . (A4) 
In order to compactify the notation, we replace each string of indices {ai, a2, • • •} by a single numerical index: 

G„,c.....a,. G[z] ; r"i"--"" := r[i] . (A5) 
We assign the choices of the variables aiQ;2 • ■ ■ a„ to the variable i using the vector 



V„ = 
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where the symbol (8) indicates the outer product. For each n, the zth component of the vector corresponds to a list 
of variables that is assigned the number i. For clarity, the results are Hsted below. To simplify the notation we drop 
the subscripts and write a Hst like "r-ir2a3" as "rra." 

[a] 2-point functions: rr ^ 1, ar ^ 2, ra — > 3, aa ^ 4 

[b] 3-point functions: rrr 1, arr —>■ 2, rar 3, aar 4, rra 5, ara 6, raa 7, aaa — > 8 

[c] 4-point functions: rrrr — > 1, arrr — > 2, rarr — > 3, aarr 4, rrar 5, arar 6, raar — > 7, aaar — > 8, rrra 9, 
arra — *■ 10, rara —^ 11, aara — )■ 12, rraa 13, araa — )■ 14, raaa — > 15, aaaa 16 



Using expressions like Ijl2ll . it is easy to show that 

ri , -prrrr---r r\ 



G[2"]=0; r[l] = 0. 



(A7) 



APPENDIX B: DISCONTINUITIES 



Each discontinuity has a pole and a cut contribution and has the general form: 

dx{K) = -2TTi ^ Zx(k) n 5{ko - nuj.j;{k)) - 2TTi(3x{ko,k) ; x ^ {T,L,+,-} . 



(Bl) 



The functions UTip), ^l{p), ^+{p) and are the roots of the denominators of the corresponding propagators. 

The expressions for these functions, and for the residues Zi and cut functions (3i are given below. 
The cut functions are: 



PT{po,r) 

PL{po,r) 

P+{ro,r) 

f3-{ro,r) 
The residues are: 



2r^m%po {pI 



^2p2 ^j.2 _ p2^2 ^ ^2j,5 _ 2p^r^ + In m^por^ + 2m'Qp'^r — In 



(B2) 



tt-^PqUi'Ij + ir^ + 2rn^r — In 



rg+r 



m%po 





2r^m'j (r — ro) (r^ 












(2 (r — ro) + m'j (r (in 


ro+r 
ro-r 


+ 2) - In 


ro+r 
ro-r 


■■»)) 



2r^m'j (ro + r) 6 (r^ — Tq) 



In 



'•0+'' 
ro— r 



- 2 



In 



rn+r 
ro-r 



ro j - 2r2 (r + Tq) j 



ZL{r) 
Z+ir) 
Z-{r) 



uj+{rY — 
2mj 



2to^ 



The positions of the poles are determined from the solutions of the equations: 

ojrir) + r 



In 



+ 2mrr - In 



WL(r) + r 



ujL{r) - r 



uiQUOLir) = 



(B3) 



(B4) 
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In 



u>+{r) + r 



w+(r) - 
2r (nij + r {r + lj_ (''))) — In 



(r — uj-^{r)) nij + 2r (rnj + r (r — a;_|-(r))) = 



uj^ (r) — r 



nij: {r + uj^{r)) = 



APPENDIX C: WARD IDENTITIES AND KMS CONDITIONS 

The ward identities at finite temperature have the same form as at zero temperature. As discussed in section ITTTI 
there are 7 components of the 3-point function and 15 components of the 4-point function, each of which has its own 
ward identity. These ward identities are expressed in a compact notation as: 



(Pout - Pin)^rhtl^[l](Pin, Pout) = Shtl[ jl ] (Pn) " Ehtl[ ] (Pout) 
Q'^MhtV[l](Pi„,Ql,Q2,Pout) = rhtl. [ Jl] (Pin, Pout -Ql) -rhtl,[j2](Pir 



(CI) 



Qi,Po, 



The index in square brackets denotes the Keldysh component of the vertex. The complete sets of indices are given in 

In the htl Umit, the ward identities have a somewhat simpler structure. We give only the results that are needed 
for this calculation (using S —P + Q): 



Q^Tm^'^[2]{Q,R)^YMV;^t{R)■ 
g^rhtl'^'^[4](Q,P) = 0; 
0^rhtl'^''[6](Q,P) = Shtl,- (P); 
g,rhtP''[2](P,Q) = Shtl^,t(P) ; 
g,rhtP''[3](P,Q) = Shti;:,„(P); 
Q,rhtP''[4](P,Q) = Ehtiryn,(P); 

Q^,Q,M^'''^[2]{Q,P,-P,Q) 

QmQ-^^''""[4](Q,P,-P0) 
QMQ-M^''"[6](g,P,-P,Q) 



p,rhtP''[2](g,p) = sht4,(Q) 

P^rhtP''[4](Q,P) 
RyTU\^"'[&]{Q,R) = 

p^rhtP''[2](p,g) = shti^,t(Q) 
p^rhtP''[3](p,g) = shti^,t(g) 
p^rhtP''[4](p,g) = 

s]hti;,,(p) + shti;,,(5) 

NF{so){YM\l,^{S)-Yht\l^,{S)) 
NF{ro){YM\l,^{R)~YM\l^,{R)) 



(C2) 



(C3) 



Terms which are contracted with the part of the photon propagator that is proportional to g^^ cannot be written 
in terms of the self-energy using ward identities. We simplify these terms using the kms conditions for the 3-point 
functions: 



rhti^[4](g,p) 



(C4) 



= -Nb (po) rhti^[ 2 ] (g, R) ~ Nf (go) rhti^[ 3 ] (g, r) + Nb (po) rhti;[ 5 ] (g, r) + Nf (go) rhti;[ 5 ] (g, p) 
rhti^[6](g,P) 

= Nf (ro) rhtl^ [2]iQ,R)-NF{qo) Thtl^ [5]iQ, R) + Nf (go) Thti; [3]{Q, R) ^ Nf {ro) Thti; [ 3 ] (g, P) 

rhti^[ 4 ](p,g) 

= Nb (po) rhtl^(2)(P, g) - Nf (ro) rhtl^(3)(P, Q) - Nb (po) rsthtl^(5)(P, Q) + Nf (ro) rsthtl^(5)(P, Q) 
Some additional relations between vertices are: 

rhti; [5]{x,Y) = rhti^ [ 2 ] (f, x) (cs) 
rhti;[2](x,r) = rhti^[5](y,x) 
rhti;[3](x,r) = rhti^[3](y,x) 

rhti^[3](p,g) = rhti^[5](g,p) 

rhti^[2](p,g) = rhti^ [3 ](g,p) = rhti^[ 2 ](g,p) 



APPENDIX D: EXPLICIT VERTEX EXPRESSIONS 



The following results are obtained by direct calculation. They are used to rewrite Eqn. I|26|) . 

Ehtl5^*(P) 



rhtl[ 2 ](g,P) 



00 



50 



(Dl) 
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rhti[2](g,i?)o. 
rhti[ 2 ](g,i?)y 

rhtl[5](i?,Q)oo 
rhtl[5](i?,Q)o. 

rhti[5](ii',Q)y 



The power counting arguments of Q indicate that for the real part of the self-energy, corrections to the result for 
the 1-loop diagram obtained by expanding to next-to-leading order in the ratio of the external momentum to the loop 
momentum, and contributions to the 2-loop diagrams from the region of the phase space that corresponds to both 
loops hard, could contribute at relative order e, when compared to the leading order contribution. These contributions 
are referred to as (1-terms) and (2-terms) in the introduction of this paper. The leading order contribution to the mass 
is of order e^T^/qg ~ eT (see Eqn. JSJ). Explicit calculation shows that contributions to the mass from (1-terms) 
and (2-terms) are not of order e{eT) = e^T but rather of order e'^lnl/e. A similar result was obtained by Schulz in 
for the gluonic plasma frequency. 

In this appendix we show why the power counting argument fails for each of the (1-terms) and (2-terms) individually, 
in our example. It is not known if power counting always over estimates the contributions of these terms. Throughout 
this appendix we use the feynman gauge. We use K to denote a hard momentum, and Q is the soft momentum (go, Q). 

Pairs of lines with the same statistics may contribute a factor Nx{k) — Nx{k -\- q) ^ q/T and thus reduce the 
contribution of a diagram by one power of e. We show below that even without this additional power of e, the 
contributions described above do not contribute at order e(eT). For the 1-loop diagram and two of the 2-loop diagrams, 
we make use of this knowledge in advance by considering only terms that depend on the bose-einstein distribution 
function. For the 2-loop bubble diagram (see Fig. 5b) it is easier to consider terms that depend on the fermi-dirac 
distribution function. Terms that depend on the bose-einstein distribution function contain a factor NB{ko){Dret{K) — 
Dadv{K)). This factor is even in ko, is proportional to the delta function S{K^), and has mass dimension -2. In a 
power counting analysis (keeping only finite temperature contributions) we can make the replacement: 

NB{kQ){Dret{K)- Dadv{K)) ^ nb{k) ' ko^nk (El) 

where the index n is summed over ±1. Similarly, for terms that contain a fermi-dirac distribution function we can 
write: 

NF{ko)iSretiK) - Sadv{K)) -^^nf{k); ko^nk (E2) 

We consider below all possible contributions from (1-terms) and (2-terms). The diagrams we need to consider are 
shown below. 



^^-'^Ehtir(i?)-. 



r'^qo \ rrij: 



^ (-5,, (^rro - ^i?2Shtir (i?) j - nf, [irro + ^ i^' ~ Sr^) S]htlS^*(i?) 
Shtlg'*"(i?) 

"i^»^-3^htir(i?)^ 



APPENDIX E: POWER COUNTING ARGUMENTS 



K .- 



Q K+Q 



Q 



FIG. 4: The 1-loop htl. 
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, K+Q 



Q K+Q K+Q Q 

(a) 



K+Q 

(b) 



Q Q 



/ Q 



(c) 



/ K 



FIG. 5: The 2-loop diagrams. 



(1) Corrections to the result for the 1-loop diagram obtained by expanding to next-to-leading order in the ratio of the 
external momentum to the loop momentum. Consider the diagram shown in Fig. 4. Power counting indicates that 
the contribution of this diagram is of order: 



diagram (1) ~ ^ dkk^ 

^—^^ phase space 



nb{k) 

fc2 



1 



if •g + Q2/2 



ko—nk 



Tr(7o7^(^ + Q hn 



ko—nk 



(E3) 



7 propagator fcrmion denominator 

Doing the trace and expanding the denominator of the fermion propagator we have: 



diagram (1) ^ 
The leading order contribution is: 



y dkkntik) -fi + l^l fl- 



n=±l 



2koqo 



numerator 



ko—nk 



diagram (l)j^ ~ dk k ribik) 
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'■II 



(E4) 



(E5) 



which agrees with ^ . The next term in the expansion gives zero when summed over n and thus the term of relative 
order e that is predicted by the power counting arguments of identically cancels. 



diagram X! dk n nb{k) — > e^T • (0) 

The largest sub-leading term is of order: 



n=±l 



diagram (1) 



nnlo 



e^qa ^ rn,{k) 



^qo dk nf,(fc) k Prin- 



e^Tln{l/e) 



(E6) 



(E7) 



where the arrow indicates the expression that would be obtained if the expansion was done carefully enough to 
preserve the term that regulates the apparent infra-red divergence in the log. 

(2) Contributions from the double-rainbow diagram with both loop momenta hard. 

Consider the diagram in Fig. 5a. Since both loop momenta are hard, we can replace the inner loop by the leading 
order piece of the 1-loop self-energy. The leading order contribution is: 



diagram (2) 



E 

n=±l 



dk k-" 

phase space 



»b(fc) 
fc2 



7 propagator 



1 



{K-Q + Qy2y 



kQ—nk 



/cq —nk 



fermion denominators 



numerator 



If we collect the powers of e it appears as if we have: an explicit factor of in front; a factor contained in the 
htl self-energy; and a factor of l/e^ from the 2 powers of 1/qo introduced by the 2 fermion propagators. By simple 
dimensional analysis this should produce a contribution e(eT) as predicted in 0. However, if we look closely at the 
numerator of this expression we see that, as in case (1) above, this leading order term cancels identically. If we do 
the trace we obtain two kinds of terms: 



ko—nk 

ill) if2shtlo(if)|fe„=„fc^O 



n k e^T^ 



(E9) 
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The first term is zero when we sum over n and the second term is zero when we substitute in for kg. 
(3) Contributions from the bubble diagram with both loop momenta hard. 

Consider the diagram in Fig. 5b. Since both loop momenta are hard, we can replace the bubble by the leading order 
piece of the 1-loop photon polarization tensor. In this case it is easier to consider the terms that depend on the 
fermi-dirac distribution function. Using (IE2|I the leading order contribution is: 



diagram (3) ~ ^ ^ ^^^^ Tr(7o7M^ 7.)5'^"n.A/\ (ElO) 

■^-^ V fc^ ^ [K ■ Q + Q'^ I ly ko=nk ko=nk 

^—^^ phase space v ^ ^ ^ ^ v 

fci 'inion denominator -y pi-opagators nuTfierator 



If we collect the powers of e it appears as if we have: an explicit factor of in front; a factor contained in the 
htl polarization tensor; and a factor of l/e^ from the 2 powers of 1/qo introduced by the 2 photon propagators. By 
dimensional analysis this should produce a contribution e(eT) as predicted in 0- However, if we look closely at the 
numerator of this expression we see that, as in cases (1) and (2) above, this leading order term cancels identically. The 
htl polarization tensor contains three different structures: gr\, KrK\/K^ and ko{groK\, g\oKr}. Doing the trace and 
contracting with each of these tensors we find that the first two produce a term proportional to fco (which gives zero 
when summed over the index n), and the third is proportional to K'^/k'^ (which gives zero when we substitute in for ka). 

(4) Contributions to the crossed-rainbow diagram with both loop momenta hard. 

Consider the diagram in Fig. 5c. Since both loop momenta are hard, we can replace the left hand loop by the leading 
order piece of the htl vertex. The leading order contribution is: 

diagram (4) ^ E ^ (K ■ qI Qy2) L.k ^^^^"^-^ ^^^^^^^^^^^l..^. ^^''^ 

^—^^ phase spaee ^ v ^ ^ ^ ^ v ^ 

7 propagator formion denominator numerator 

There are not enough powers of 1/qo to produce a contribution of order e(eT). 

The conclusion is that, for the case of a static fermion, the power counting argument of 121] over estimates the 
contribution of: (1) the subleading corrections to the 1-loop htl result and; (2) each of the 2-loop diagrams with both 
loop momenta hard. In every case, the nlo terms predicted by are present, but they cancel exactly at the relative 
e order that is being considered. 
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